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PREFACE 

This  Is  the  first  of  a  set  of  three  reports  on  the 
theoretical  analysis  of  the  shock  wave  produced  by  an  tinderwater 
atomic  explosion.   The  purpose  of  the  analysis  is  to  determine 
the  peak  pressure,  Impulse  and  arrival  time  of  the  shock  as  a  func- 
tion of  distance  from  the  explosion  and  of  the  energy  released.   In 
addition  the  pressure  and  velocity  distributions  behind  the  shock 
will  be  found. 

The  various  quantities  mentioned  above  have  been  deter- 
mined in  two  different  ways.   These  are 

1.  Niimerical  solution  of  the  equations  of  compressible  fluid  flow, 
carried  out  on  the  UI^IIVAC  at  New  York  University, 

2,  Approximate  analytical  solution  of  the  equations  of  compressible 
fluid  flow. 

The  analytical  solution  is  partially  described  in  the 
present,  report,  ind  the  description  will  be  completed  in  the  next 
report.   The  numerical  solution  will  be  described  in  the  third 
report  of  this  set  where  its  results  will  be  compared  with  the 
analytical  results.   It  will  be  shown  that  the  two  different  methods 
of  solution  yield  results  in  good  agreement  with  each  other.   Thus 
the  possibility  of  solving  this  problem  by  an  analytical  method  will 
be  demonstrated.   Therefore  it  will  be  possible  to  solve  additional 
problems  of  this  kind  entirely  analytically.   This  will  afford  a 
considerable  saving  in  time  and  money  over  the  numerical  method,  and 
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provlde  a  better  understanding  of  the  physical  phenomena  involved. 

The  first  step  in  both  the  numerical  and  analytic 
methods  is  the  calculation  of  the  "strong  shock"  solution.   This 
solution  desceibes  the  Initial  behavior  of  the  shock.   Although 
it  has  been  studied  before,  It  is  discussed  anew  in  the  present 
report  because  of  its  basic  role  in  both  methods.   Furthermore 
various  procedures  are  developed  for  simplifying  the  calculation 
of  this  solution.   First,  of  course,  the  problem  is  formulated 
and  the  equation  of  state  of  water  is  discussed.   Afterwards  the 
equations  are  converted  to  dimensionless  form. 

The  table  of  contents  indicates  what  the  next  two 
reports  will  contain. 
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SPHERICAL  V/AVE  PROPAGATION  IN  WATER 

§1.   The  Equation  of  State  of  ^Jater 

V/hen  the  pressure  is  less  than  fifty  kilobars,  a  moder- 
ately good  fit  to  the  experimental  pros sure -volume -temperature 
values  Is  given  hj   the  relation  .  ~  ■ 

(1.1)  (F  +  P)V  =   BT   . 

In  order  to  get  the  best  agreement  with  experiment  it  is  desirable 
to  make  the  quantities  P  and  B  depend  on  the  temperature,  but  for 
the  purposes  of  this  discussion  vie    shall  assume  they  are  constant. 
The  quantity  P  is  called  the  internal  pressure  and  is  approximately 
equal  to  3000  atmospheres. 

The  enthalpy  of  the  water  is  defined  by  the  relation 


(1.2)  i  =  e  +  pV 

where   e    is   the   internal   energy.      It    cmi  be    shovxn  from  the   laws   of 
thermodynarnlcs   that   the   change   in  enthalpy  between  two   states 
(p    ,T    )    and    (p,T)    can  be    calculated  from  the   relation 

(1.3)  i    -   Iq  =      /    Cp(T',p^)dT'    +    I    iv(p,T)    -   T  m§l^|dp> 

T  pi  ^ 

where  c  (T',p  )  is  the  specific  heat  at  constant  pressure  p  . 


•*  -.  C3  .  ''      Av  ri'-*       -> 
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It  is  easily  shov/n  that,  vmen  the  state  equation  has  the  form  (1.1), 
the  second  term  vanishes  identically  and  we  can  write 


(l.U)  i  -  io   =   ^p(^  -  ^^^o^ 


where  c   is  the  average  value  of  c  (T',p  )  in  the  ranee  T  <  T'  <  T 
p  °  p    '^o  &   o  —    — 

An  equivalent  form  of  (l.ii)  is 


(1.5)  i  -  i   =  -^ 


j  (p  +  p)v  -  (p^  +  p)v^  r 


where  we  have  introduced  the  ratio  of  specific  heats  at  constant 

prossui'^e  and  constant  volume,  y  =   c   v/c   ,  and  the  thermodynamic 

fact  that  B  =  c   -  c   .   For  xi/ater  y  has  an  average  value  of  7«l5 
p    V  '  ^  '   -^ 

and  c..  is  approximately  equal  to  3«96iilj.  joule s/gm.deg. 

The  expression  for  the  internal  energy  is  easily  computed: 


(1.6)  ^'  -  ^o  "  7^1  [(p+P)v  -  (Po+P)V^]  +   P(V-V^) 


We  also  compute,  from  the  first  law  of  thermodynamics, 
the  entropy  change 

T  p 

(1.7)  s-s^  -  I   o^^L  .  J    v(p',T)^  =  Cp  m  f-  -  B  m  £i| 


p.+p 


T  p 

o  '■  o 


F-rom  this  expression  and  (1.1)  v/o  derive  the  cntropic  form  of  the 
state  equation: 
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s-s 


(1.8)  (p+P)vY  =   (Po+P)vJ  oxp( ^) 

The-  sound  speed  c  in  water  is  easily  obtained  from  its 
definition  by  use  of  (1.8),   It  is  given  by 


(1.9)  c^   =   -  v2(^)3  =  y(p+P)V  =  yBT 


§2.        The  Equations  of  Fluid  Motion 

It  is  convenient  in  discussing  spherical  blast  waves  to 
introduce  a  Lagrangean  coordinate  h  which  is  proportional  to  the 
mass  contained  inside  a  sphere  of  radius  r,   VJe  define 


r(h,t) 
(2.1)  ^     ^      I  V  ^^'^^ 


which  in  differential  forra,  using  subscript  notation  for  differentia- 
tion. 


(2.2)  V  =  r^v^ 


expresses  the  law  of  mass  conservation. 

Newton's  law  of  momontuifi  conservation  can  then  be  written 


(2.3)  u^   =   -  r\ 


VJe  shall,  moreover,  assume  that  all  motions  behind  a  hvdro- 
dynamic  shock  take  place  adiabatically,  so  that  the  law  of  energy 
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conservation  can  be  written  in  the  equivalent  forms 


(2.i+a)  S^   =  0 


or 


(2.i+b)  Vp^  +  y(p  +  P)V^  =  0 

To  make  this  system  of  partial  differential  equations  com- 
plete we  must  add  the  kinematic  equation 


(2.5)  ^t   "  "^ 


and  one  of  the  state  equations  (1.1)  or  (1.8)  . 

§3.   The  Shock  Conditions 

If  a  shock  or  pressure  discontinuity  e:cists  in  the  water 
medium,  then,  as  a  consequence  of  the  conservation  laws,  the  states 
on  the  two  sides  of  the  shock  are  connected  by  the  Rankine-Hugoniot 
conditions,  ^^7hich  may  be  ^^Jritten 


(3.1)  V^(U  -  u^)   =  V-^U 


(3.2)  V  (p^  -  Pq)   =  u^U 


r^fivc 


viJ} 
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{3»3)  ^1  ■  ^0   "   (1/2)(Vq  +  V^)(Pi  -  Pq)   • 

In  these  equations  the  subscript  1  refers  to  the  state  just  behind 
the  shock  and  the  subscript  0  refers  to  the  ambient  or  rest  state 
ahead  of  the  shock.   V/e  assvime  that  the  water  is  at  rest  ahead  of 
the  shock.   The  quantity  U  =  U(t)  is  the  shock  speed  and  we  have 
the  relation 

(34)         ■    U  =  dR/dt 

where  r  =  R(t)  is  the  position  of  the  shock  at  time  t. 

It  is  convenient  for  our  purposes  to  transform  the  equa- 
tions (3»1)  to  (3»3)  by  use  of  the  results  of  §1.   In  addition  vie 
define  the  shock  pressure  ratio 


p,+P 
^^•^^  ^1   =  ^ 


PO+P 


Then  the  shock  conditions  can  be  put  in  the  form 


(3.6)  V-L  =  V^iii^t,-^  +  l){%n,  +   P-^)"^ 


(3.7)  u^^  =  U(l  -  ^^)(^;l  -  l)(-^.i  +  t^^)"^ 


(3.8)  P-L  +  P  =  U^(l  -  ti^)^^(£,^  +  [i^)'\^ 


—  r 
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where  \i     =  (y-1)/(Y"''1)  • 

In  particular  ue    are  interested  in  the  limiting  forms 
of  (3.5)  -  (3»8)  when  the  shock  strength  becomes  very  large.   The 
resulting  relations  are  known  as  the  strong  shock  conditions: 


(3.9)  V3_  =  ia\ 


(3.10)  u^  =   (1  -  u^)U 


(3.11)  Pi  +  P  =  (1  -  i^^Ju^A 


itj.,    Taylor  VJaves 

An  important  type  of  solution  of  the  hydrodynamic  equa- 
tions and  boundary  conditions  has  been  discussed  by  G.  I.  Taylor 
[1].   The  solutions  are  also  variously  called  progressing  waves, 
strong  shock  solutions,  point  blast  solutions,  or  similarity 
solutions,  cf.  [2,3,L|.].   They  may  be  described  in  the  following  way. 
Let  r  =  R(t)  denote  the  shock  path  and  U(t)  =  dR/dt  the  shock  ve- 
locity.  We  then  ask  Xifhether  the  Euler  radius  of  a  particle  r(h,t), 
the  particle  velocity  u(h, t),  the  specific  volume  V(h,t)  and  the 
pressure  p(h, t)  can  be  represented  as  a  product  of  a  function  of 
the  time  by  a  function  of  a  similarity  variable  x,  which  depends 
on  both  h  and  t  in  some  prescribed  l^ray.   For  our  purposes  it  is 


■:  l-i-       •  l._       i- 


^v^V't''-      ''J^V     ^r::\-:<:":,Q    x\ 


';    :>i;v 
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convenient  to  raake  the  choice 


(i|..l)  r(h,t)   =  R(t)f(x,Y) 


(1|.2)  u(h,t)   =   (l-ii^)U(t)g(z,Y) 


{k'3)  V(h,t)   =  tx^V^T(x,Y) 


{k*k)  p(h,t)  +  P  =  (1  -  ^l^)-^V(t)tT(x,Y) 


where  x  =  h/H(t)  and  H(t)  =  R^(t)/3V^   . 

In  order  to  understand  the  restrictions  that  are  implicit 
in  this  form  of  solution  vie  shall  exi"uviine  the  energy  balance  in  the 
blast  wave.   Since  the  motions  behind  the  shock  are  assvuned  to  take 
place  adiabatically,  the  total  energy,  internal  plus  kinetic,  of 
the  fluid  inside  the  shock  sphere  at  any  instant  must  equal  the 
yield  energy  Q  put  into  the  blast  plus  the  internal  energy  of  the 
arabient  fluid  taken  into  the  shock  sphere,  i.e. 

1 


(/^.5) 


0„  =  it.TT  H(t)   Me  -  e^  +  I  n.^)dx 


But,  ^^rith  our  assumptions  about  the  solution. 


-Syr'     1     >    ■■ 


.-■t  »■ 
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2  2 
e-e.  +  i  u^  =  ^^'^    ^      u2(t)[Tr(x,Y)T:(x,Y)+  g^(x,Y)] 


'o   2  "       2 


2. 


+  PV  [ix'-t(x,y)-1]  -  ^•'-"^  ^  (Pq  +  P)Vq 

2p, 


and  hence 


Q  =  2TT(l-ix^)^H(t)u2(t)I(Y)  -  2Tr^l"^(l-^L^)(p^^■P)V^H(t) 


where  we  have  set 

1 
(I|..6)  I(y)  =  f   [tt(x,y)t(x,y)  +  g^(x,Y)]dx 

and  have  made  use  of  the  fact,  which  follows  from  (2.2),  that 


X 

J      [^^'r(x,Y)  -  l]dx  =  0   . 


If  life  modify  the  form  of  the  last  term  by  using  the  relations  {3*5) 
and  {3'^)}    the  expression  for  Q  becomes 

(U.7)    Q  =  2Ti-(l-tJL^)^I(Y)H(t)u2(t)[l-lAi2l(Y)('5-L(t)+p.^)] 

Since  the  yield  Q  is  independent  of  the  time,  the  expression  on  the 
right  hand  side  must  be  also.   This  can  certainly  be  arranged  by 
setting 


-■<; 


)  I  ■  .'■ 


)•  :i 


CIl       '! 


{  1  j-ii'./)' 


-■J:i  ■.  ,  , 


■■.R 
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(1^.8)  U^(t)R^(t)   =   A^ 


where  A  is  a  constant  and  by  taking  ^-i(t)  so  large  that  the  second 
term  is  essentially  zero.   The  first  requireinent  implies  that 


(1]..9)  R^(t)   =  4?  A^t^ 


and  the  second  implies  that  our  solution  will  only  be  valid  when 
the  shock  pressure  ratio  is  very  large.   The  appropriate  boundary 
conditions  for  the  Taylor  waves  are  therefore  the  strong  shock  re- 
lations (3.9)  -  (3.11). 

Wlaen  the  last  term  in  the  bracket  of  equation  (li.?)  is 
set  equal  to  zero,  the  constant  A  can  be  determined  in  terms  of  the 
yield  Q,  the  ambient  specific  volume  V   ,  and  the  adiabatic  exponent 
Y  by  means  of  the  relation 

(14-. 10)  A^   =   3V^q/2tt(1  -  y.^)h{r)       . 


0.?     < 
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s5«    Scalinp;  of  the  Equations  to  Dimenslonless  Form 

We  can,  of  course,  measure  the  physical  quantities  in- 
volved in  our  results  so  far  in  any  consistent  dimension  system. 
There  is  a  choice  of  units,  however,  x-;hich  is  related  to  the  para- 
meters of  the  problem  and  which  leads  to  results  of  much  greater 
mathematical  simplicity.   This  possibility  depends  on  the  fact  that 
all  explosion  phenomena  are  similar,  if  the  motion  is  adiabatic 
behind  the  shock,  by  which  we  mean  that  they  depend  only  on  the  ini- 
tial energy  or  yield  of  the  blast  Q,  and  upon  the  quantities  charac- 
terizing the  ambient  state  ahead  of  the  shock,  for  example  the 

pressure   p      and   specific   voliime   V      .      Vie    shall   choose    our  dimension 
•^o  -  o 

units  so  as  to  make  our  results  independent  of  these  quantities. 

They  are  as  follows: 

Length  unit:  \,    where  \-^  =  Q/p  +P 

Velocity  unit:   c  ,  vjhere  c   =  y(p  +P)V 
■^         o         o      o    o 

Time  unit :  x/c 
o 

Pressure  -onit:   p   +  p 

Energy  unit:   c'^ 

Mass  unit:   X^/v 

o 

Specific  voliame  unit:   V 
'■  o 

Our  previous  results  can  then  be  restated  in  the  form 


(5.1)  1  =  -~  i^V  (enthalpy) 


i-/here    "^  =   (p   +  P)/(p     +  P)    is   the   pressure   ratio. 


:tX   h. 


,x.-.;    ;c     , 


::   ,k:- ':'-■■: 


y\Ui:i     1...   .-:;■ 


'  '   "■  ■-•;".T       ■'" 


',  ..'        v<^      U.L  ;:  .!- Y 


■".  .,  -i       "i  1 
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IS.2) 


Y~  ( — r-   t^    +  k)V       (internal  energy) 


where  k  =  P/(p   +  P), 


(5.3) 


^V^  =  exp  3., 


(entroplc  state  equation) 


where  5.,  =  (S-3^)/c    * 


(5.U) 


=  t,V  =  B,/ 


(sound  speed) 


where  B,..  =  y(y  -  l)c  /c 


(5.5) 


V  =  r~r. 


h 


(continuity  equation) 


(5.6) 


^t  =  ^^ 


(Lagrangean  velocity  definition) 


(5.7) 


„  -    -12, 


(momentuin  equation) 


(5.8)       ^^t  "*"  '''^^^-  ~  *^       (adio.batic  energy  equation) 


(5.9)       V^(   ^  +  p.^)  =  l.L^  ^  +  1    (first  shock  condition) 


(5.10)      ^l^^l  "*"  ^   ^    =  U(l  -  n,^)(^,  -  1)   (second  shock  conditi 


on. 


2^..2 


(5.11)       ^-j_  +  tx^  =  (1  +  tx^)U 


(third  shock  condition) 


C  J. 


I-.:-' 


1     '       /  - 


H    * 


J 
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p 

(5»12)     v.,   =  p.  (first  strong  shock  condition) 


(5»13)     u,   =   (1  -  ii"^)!!       (second  strong  shock  condition) 
(5»ll4-)     v>-|   ~   (1  "^  I-*-  )U'"      (third  strong  shock  condition)  . 

In  order  to  get  a  feeling  for  the  magnitude  of  the  quan- 
tities involved  we  note  that  if  vje  take  p   +  P  to  be  3000  bars,  the 
unit  of  length  X  will  then  be  10  meters  if  the  yield  is  'J\*o   r.ietric 
tons  of  TNT  and  x-rill  be  100  neters  if  the  yield  is  71*6  r.ietric  kilo- 
tons  of  TNT  or  its  equivalent.   If  the  yield  is  20  kilotons  of  TNT 
equivalent,  the  length  unit  is  65.3  ■-■leters. 

16.    The  Solution  of  R.  Primakoff  for  y  =  7. 


If  the  equation  of  state  has  the  form  (1.8)  with  Y  =  7 
and  [J,  =  3/i|  an  explicit  solution  of  the  dynra;!ical  equations  (5*5) 
-  (5*8)  can  be  given  cf.  [2,ii.],  which  is  co.'npatible  x-xith  the  strong 
shock  conditions  (5 'IS)  -  (5«lU)»   (This  special  solution  was  also 
discovered  earlier  by  VanVleck  but  its  significance  for  water  blasts 
v/as  not  noticed.)   It  is  a  special  case  of  the  progressing  wave  or 
point  blast  solution  discussed  by  G.  I.  Taylor.   As  q   preliminary'- 
we  compute  the  total  mass  of  vmter  H  which  is  contained  inside  the 
shock  sphere  at  time  t.   Aside  from  a  factor  [[.Tr  it  is  given  by 

R(t)  R(t) 

H(t)   =  {  V"^r^dr  =   /     ^o'^^o 
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where  r  is  the  initial  distemce  of  a  particle  from  the  center  and 
the  second  relation  folloi%rs  frora  the  equation  of  mass  conservation, 
But  in  the  latter  form  the  integration  can  be  perforraed  explicitly 
with  the  result 

(6.1)  H(t)   =  i-V3 

This  equation  is  independent  of  the  value  of  y  • 

The  Primakoff  solution  for  y  ~  7  can  be  written  in  the 
form 

(6.2)  ?(h,t)   =  R(t)x°-^^ 

(6.3)  u(h,t)   =  0.2i;u(t)5^'2^ 
(6.L)  V(h,t)   =  0.75:"  '^''^^ 

(6.5)  -?(h,t)   =  7Ay'"(t)p/^  1.75U^(t);"°-7^ 
where 

(6.6)  X  =  h/n(t) 

The  values  of  these  quantities  just  behind  the  shock  are  foi:^d  by 

setting  X  =  1. 
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It  is  also  instructive  to  v/rite  these  quantities  as  func- 
tions of  the  Euler  distance  r.   VJe  define 

(6.7)  y  =  ?/R(t) 

and  then  we  have 


(6.8)  h(r,t)   =  H(t)y^ 


(6.9)  u(r,t)  =■■     0.25  U(t)y 


(6.10)  v(?,t)  =  0.75  y  "■"■ 


(6.11)  t^?>t)  =  1.75  U^y^ 

These  relationships  are  shown  in  Figure  1. 

The  shock  speed  and  radius  are  determined  since 
I(y=7)  =  ^4-/3  by  application  of  ()4..6),  from  the  relations 

(6.12)  U(t)   =   (18/7tt)°*%  "^-^ 

(6.13)  rut)   =   (225/ll|.'(T)^'^t°*^ 

The  folloxiing  formulae,  which  we  make  use  of  later,  follov;  from  (6.1), 
(6.12)  and  (6.13): 


v-  .  ■  ; 
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(6.110  H  "^dH/dt   =  3R  U  "-^   or  dS/dt   =  R^IJ 


(6.15)  dU/dt  =  1.5  R  ''^^^ 


The  use  of  the  results  derived  in  this  section  may  be 
illustrated  by  the  following  exaraple.   Suppose  \-jq   know  that  for  a 
certain  explosion  in  water  the  total  energy  going  into  the  blast 
wave  was  the  equivalent  of  10  metric  Icllotons  of  TKT.   VJe  wish  to 
find  the  conditions  existing  at  the  time  the  shock  sphere  has  a 
radius  of  20  meters. 

VJe  may  assurae  without  making  any  significant  error  that 
v;ater  has  the  value  of  y  ~  7  and  p  +  p  =  300O  bars.  In  order  to 
find  the  unit  of  distance  \-ie   use  the  relation 

^  ~  yO/Po^  =   /l0xl+.  185x10-^ '^dyns-ciV3.000xl0'^d3^e/cm^ 


=  Vl-395xl0"^^cra^  =  Vl39.5xl0^  meters  =•  51.8  meters 

Hence,  the  dimensionless  shock  radius  is  R  =  20/51.8  =  O.386  .   The 
so\ind  speed  in  the  ambient  region  is 

=  /yV~Tp~+pT  =•  v/7xlcm-^/gmx3xlo'^ dynes/cm^  =  /2.1  x  10^  m/s 


=  II4.50  m/sec 


'  :-;■.'■.    ' 
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Prom  equation  (6.12)  we  compute  the  dimensionless  shock 
speed,  U  =  /i8/7ttR^  =  /l8/7Tr  ( O.382  }^  =  3.79  .   The  actual  shock  soeed 
is  3.79  X  lq.50  =  5^00  Ill/sec. 

Prom  equation  (6.11)  v;ith  y  =  1,  vje  can  compute  the  dimen- 
sionless  shock  pressure  ratio,  ^-  =  1.7p  U  =  25«2  .   (if  this  value 
had  been  less  than  10,  the  use  of  these  simplified  formulas  would 
not  be  justified.)   The  actual  pressure  at  the  shock  is 
p  =  25.2  X  3000  -  3OCO  bars  or  p  =  72.6  kilobars,  assuming  that  p 
is  approximately  zero.   The  pressure  at  a  point  10  meters  from  the 
center  of  the  shock  sphere  is  fo-ond  from  (6.7)  and  (6,11): 
y  =  r/R  =  r/R  =  10/20  =  0.5,  and  ^(0.3', t)  =  25.2  x  (0.5)^  =  3.13' 
or  p  =  6,1].5  kilobars. 

It  is  perhaps  appropriate  at  this  point  to  say  something 
about  the  range  of  applicability  of  the  results  obtained  so  far. 
Two  principal  simplifying  assumptions  have  been  made,  one  hydro- 
dynamic  in  nature  and  the  other  theiTiiodynamic . 

The  C'lief  feature  of  the  hydrodynai'iic  sim.plification  vjas 
the  restriction  to  strong  shocks.   Thus  we  neglected  consistently  term 
of  order  \/'%-.     .   A  crude  estimate  would  tharefore  put  the  error  due 
to  this  assumption  in  the  sai"ie  order  of  magnitude.   If,  for  example, 
a  maxlmxim  error  of  ten  percent  is  allov;able  vje  sbould  not  use  the 
Primakoff  solution  when  the  pressure  behind  the  shock  has  droprjed 
below  about  3O5OOO  atmospheres. 

The  errors  due  to  the  thermod^ynamic  assuraption  that  the 
equation  has  the  simple  Tait  form  (1.8)  are  much  more  difficult  to 
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gauge,  since  our  knov/ledge  of  the  t he rrao dynamic  properties  of  water 
is  very  meager  for  high  pressures  and  temperatures.   In  one  respect, 
however,  it  is  possible  to  estimate  the  effect  of  a  change  in  the 
equation  of  state.   It  will  be  recalled  that  the  Primakoff  solution 
holds  only  for  the  value  Y  ~  7>  vjhile  the  average  value  of  y  fo^  water 
was  stated  earlier  to  be  about  7 '15  • 

In  the  succeeding  sections  v;e  will  sIlovj  hovi   the  Primakoff 
solution  can  be  modified  by  perturbation  procedures  to  give  the  solu- 
tion when  the  shock  strength  is  not  large  and  when  y  is  different 
from  the  value  seven. 

g 

s? .    The  Small  Y"7  Approximation 

In  order  to  improve  the  Primakoff  solution  of  the  previous 
section,  so  that  it  holds  for  values  of  y  1^"^  the  neighborhood  of  7, 
we  return  to  the  general  equations  of  fluid  motion  (3.5)  -  (5«8)  and 
the  corresponding  strong  shock  relations  (5.12)  -  (5.1i|-).   l^e  replace 
Y,  whenever  it  occurs  in  the  equations,  by  7  +  e  and  make  the  assump- 
tion that  for  small  enough  values  of  e  the  quantities  r,  u,  V  and  p 
have  a  convergent  power  series  expansion  in  e  of  which  the  first  term 
is  the  Primakoff  solution.  Moreover,  we  look  for  progressing  wave 
or  similarity  type  solutions  in  which  the  higher  order  terms  in  the 
e-expansion  depend  on  h  and  t  in  the  same  way  as  in  the  Primakoff 
solution,  i.e.  they  can  be  written  as  the  product  of  a  function  of 
the  time  by  a  function  of  x,  where  x  =  h/H  and  H  denotes  the  mass  of 
the  perturbed  shock  sphere.   More  precisely,  we  set 
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r(h,t)  =  ?  +  eS  x^-^^f  (x)  =  K(t)x°'^^  I  (R/K)°*'^^  +  ef(x)  j 

u(h,t)  =  a  +  0.25etJx°*2^g(x)  =  0.2SU(t)x°-^^  Ue/^?'^^   +  eg(x)| 
(7  1) 

V(h,t)  =  V  +  0.75ex"°'^^T(x)  =  0.75"°*^^  |(R/H)^'2^  +  eTCx)  } 

fC(h,t)  =  1  +  l»75etj2x°*'^^7i(x)  =  1.75u^(t)x°*'^^/(R/R)^-^^  +  eTt(x) 

where  we  have  retained  only  two  terms  in  the  e-series  expansion. 
But,  since  R(t)  =  r(H,t)  and  U  =  dR/dt,  we  have 

(7.2)     R/R  =  (R/K)^*'^^  +  ef(l)  =  1  +  l4-ef(l)  =  U/U   . 

I 

Consequently,  relations  (7.1)  may  be  written 

.25  r, 

a.k)  u(x,t)  =  0.25U(t)x°'^5  K  ^   e[g(x)  -  f(l))i 

(7.5)  V(x,t)  =  0.75x"°"2^  .fl  +  e[T(x)  -  3f(l)]l 

(7.6)  ^(x,t)  =  1.75U^(t)x°'''^  \l   +  e[7i(x)  +  f(l)]j 

The  values  of  these  quantities  just  behind  the  shock  may 
be  obtained  by  setting  x  =  1  in  (7.i|-)  -  (7*6),  so  that 


(7.3)     r(x,t)  =  R(t)x°'''^  il  +  e[f(x)  -  f(l)] 
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(7.7)  u^(t)  =  0.25U(t)  jl  +  e[g(l)  -  f(l)](  =  (1  -  li^)U(t) 

(7.8)  V^(t)  =  0.75  jl  +  e[r(l)  -  3f(l)]  I  =/ 


•  r^  n^i4-\       1  J.  ^f  ttM  ^  J.  f  M  \  1  f  =  -"■"     2\^,2, 


(7.9)    q(t)  =  1.75  U^(t)   1  +  tlMl)   +  f(l)]f  =  7(1  -  [I   )li   it)    , 


where  the  last  equality  is  just  the  statement  of  the  strong  shock 
relations.   Since,  by  a  simple  calculation. 


(7.10)   M-^  =  0.75(1  +  e/2k),    (1  -  ^i^)  =  0.25(1  -  e/8) 


we  get,  by  comparison,  the  boundary  relations 

(7.11)  til)   -  3f(l)  =  IM 

g(l)  -  f(l)   =  -  1/8 
7l(l)  +  f(l)   =  -  1/8 

Similarly  the  insertion  of  (7.3)  -  (7.6)  into  the  equa- 
tions of  fluid  motion  (5.5)  -  (5.8)  gives  the  differential  equations 
satisfied  by  the  perturbation  functions 

(7.12)  t(x)  -  3f(x)   =  l+xf'(x) 

(7.13)  g(x)  -  f(x)  =  -  12xf(x) 

(7.11+)   3ti(x)  +  i|XTi'(x)  =  3g(x)  +  l+xg'(x)  -  6f(x) 


i^^ 


{ .t )  'j' '  "ii  *■ 


n 


jKii'J    -    (.OaJ3    +   if   (t]'V'lS,0    ^-    (:^    j;tf 


I,  >. 


■■'    ( i-  ^  1  j  3 


t  a    r 


1  i  ,i)-'.  ^  (i;:i 


'■•C       f 


.  .''Z.i  •=.     ': 


t'-    ■'i^i;\-:.^7.^.^t    .?;f;'     i'r.;;:     ai    v 


X]  '.. 


:i.;i    v'~:j'L;'3i'0iif    i>il:.     ^,ri 


fl^T? 


V     * 


anu:'  .:'•?. :P  ^   a;  r  jn^' 


/  r 


L-' 


^      (x)-i:   *    .:0t 


(■•0*"ix>x  ~   -     (xi""*      -    !.'x)§ 


.  \    rf. 


\K>     .  /.  ,.     4 


t\-)^L 


-22- 

(7.15)  i4.x7i'(x)   =   1  -  28xr'(x) 

The  last  of  these  equations  can  be  integrated  immediately 
to  yield 


(7.16)         ii(x)  +  Itix)      =  ^  In  X  +  20  f(l)  +  i 


which  combines  with  the  other  equations  to  give  a  single  second  order 
equation  for  the  function  f(x): 

(7.17)  I6x^f"(x)  +  50xf'(x)  +  l5f(x)  =  1^  In  X  +  l5f(l)  +  | 
whose  solution,  subject  to  the  boundary  conditions  (7»11),  is 

(7.18)  f(x)  -  f(l)  =  1^  m  X  +  ^  x-^/^  -  ^ 
From  this  solution  we  obtain  very  easily 

(7.19)  r(x)  -  3f(l)  =  ^  In  X  +  ^  x-^/^  +  ^ 

(7.20)  g(x)  =  f(l)  =  1-  m  X  ^  ^1  x-^/S  .  ^ 

(7.21)  Tt(x)  +  f(l)   =  -  ^  In  X  -  ^  x'^^^  -  ^ 


The  solution  is  completed  by  computing  the  relation  between 
U  and  R,   This  is  done  most  easily  by  use  of  the  results  of  sij..   In 
particular  we  use  equations  ([{..8),  (l4.,10)  and  (U.6)  to  obtain  the 
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result  that 


(7.22)  uV  =  A^  =  i|  (1  +  eA)t 


from  which  It  follows  that 

(7.23)  R^-^(t)   =  -i^  (1  +  e/8)t 

A  further  by  product  is  the  evaluation  f(l)  =  l/60,  although  this 
result  is  only  useful  for  a  comparison  of  the  Taylor  solutions  for 
Y  =  7  and  y  =  7  +  e  . 

The  final  expressions  for  r,  u,  V  and  e  as  functions  of 
the  variables  x  and  t  can  be  written 

(7.2i|)     r/R(t)  =  x°-2^[l  +  |-(1  m  X  +  i^  x-°-^25  .  1^^^ 

(7.25)  u/U(t)  =  \   x°-25[i  ,  |_(1  1,  ,  ,  17  ,^0.625  .  k|)j 

(7.26)  V  =  \  x-0-2%1  .  |_(3  1^  ^  ^  1_  ^-0.625  ^  ^j  j 

(7.27)  C/U^  =  I  x°-7^[l  -  1^(1  m  X  f  7_  ^-0.625  ^  3|) ^   . 

For  the  understanding  and  application  of  these  results 
it  is,  however,  more  convenient  to  express  the  results  as  f\inctions 
of  the  Euler  radius  r  and  the  time   .   The  inversion  of  (7.2i|)  - 
(7.27)  in  terms  of  the  variable  y  =  r/R(t)  gives 
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(7.28)    X  =  h/H(t)  =  y^[l  -  Idn  7  +  j^   y"^*^  -  J^)  ] 


(7.29)  u/u^(t)  =  y[l  +  ^(y"^*^  -  1)] 

(7.30)  V/V^     =  y"^[l  +  ^ik   m  y  +  Jq  j'^'^   -  Jq)] 


(7.31)       C/4i(t)  =  y^[l  -  1^(1+  m  y  +  i^  y"^*^  -  ^)  ] 


The  values  at  the  shock  are  given  in  (7.7)  -  (7.10). 

The  terms  multiplied  by  e  in  the  relations  (7.28)  -  (7.31) 
are  plotted  in  Figure  2.  The  corrections  are  probably  valid  even 
for  moderately  large  e  within  the  outer  half  of  the  shock  sphere 
but  should  not  be  taken  too  seriously  in  the  neighborhood  of  the 
center. 

By  using  the  last  of  equations  (7.1)  It  is  very  simple 
to  calculate  the  effect  of  a  change  of  y  o^  the  shock  pressure 
curve .   Thus 


^l(t;Y)/  li(t;7)  =  (R/R)^-^^  +  eit(l) 


=  1  +  (5f(l)  +  niljix   -  7) 

=  1  +  (io  -  ^)(Y  -  7)  =  1  +^^(7  -  Y) 


Hence,  we  have,  for  example, 
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?-^(t;5)   =    1.15  fi(t;7) 
^^(t;6)   =   1.075  li(t;7) 
f;^(t;8)   =   0.925  li(t;7) 
£^j_(t;7.l5)   =   1.013  li(t;7) 

It  should  be  emphasized  that  the  results  above  apply  to 
equal  times  on  the  shock  pressure  curve.   If  the  comparison  is  made 
at  the  same  shock  radii,  the  shock  pressure  ratio  ^,  for  y  near  7 
differs  by  a  factor  (1  +  •^-)  from  the  Primakoff  pressure  ratio, 
i.e.,  it  is  greater  when  y  >  1   and  less  when  y  <   1.      This  is  just 
the  opposite  of  the  results  shown  above  and  is  explained  by  the 
fact  that  the  shock  radius  versus  time  curve  is  affected  by  the 
change  in  y*  ^^  is  shown  by  formula  (7.2),  which  can  be  written 

R(t;Y)/  ^(t;7)   =   1  +  ^ 
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